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Abstract 

We consider a scenario where a supersymmetric model has multiple dark matter particles. Adding 
a U{1)' gauge symmetry is a well-motivated extension of the Minimal Supersymmetric Standard 
Model (MSSM). It can cure the problems of the MSSM such as the //-problem or the proton decay 
problem with high-dimensional lepton number and baryon number violating operators which R- 
parity allows. An extra parity ([/-parity) may arise as a residual discrete symmetry after U{1)' 
gauge symmetry is spontaneously broken. The Lightest [/-parity Particle (LUP) is stable under 
the new parity becoming a new dark matter candidate. Up to three massive particles can be stable 
in the presence of the i?-parity and the [/-parity. We numerically illustrate that multiple stable 
particles in our model can satisfy both constraints from the relic density and the direct detection, 
thus providing a specific scenario where a supersymmetric model has well-motivated multiple dark 
matters consistent with experimental constraints. The scenario provides new possibilities in the 
present and upcoming dark matter searches in the direct detection and collider experiments. 



I. INTRODUCTION 

Current cosmological data [l| indicate that about 23% of the energy in the Universe is 
in the form of cold dark matter (CDM). The origin and the nature of the CDM is one 
of the biggest puzzles in theoretical particle physics today. Since all known particles are 
ruled out as dark matter candidates, the dark matter problem provides one of the strongest 
phenomenological motivations for new physics beyond the Standard Model (SM). Among 
the many possibilities, supersymmetry (SUSY) is perhaps one of the best motivated new 
physics scenarios: it resolves the fine-tuning problem of the SM, and may provide natural 
dark matter candidates among the spectrum of new particles. 

It is not SUSY itself that guarantees the stability of the lightest supersymmetric particle 
(LSP), but rather an additional Z-i symmetry (i?-parity ^) which is introduced to solve the 
proton decay problem: without i?-parity or some similar type of stabilization mechanism, 
the superpartners of the SM fermions would be expected to mediate proton decay at an 
unacceptably high rate. -R-parity has attractive features. First, it protects the proton from 
decaying via renormalizable lepton number violating operators and baryon number violating 
operators. Second, it allows TeV scale SUSY to evade the stringent bounds from electroweak 
precision data. Finally, in the presence of i?-parity, the LSP is stable and can be a viable 
dark matter candidate. For all these reasons, the Minimal Supersymmetric Standard Model 
(MSSM) with i?-parity has been the most popular supersymmetric extension of the SM, and 
most SUSY dark matter studies have been confined to this setup jsl. Within the MSSM, the 
lightest neutralino (a mixture of neutral gauginos and higgsinos) has been the only viable 
CDM candidate since the other possibility (the left-handed sneutrinos) has been ruled out 
by combining two constraints for a viable dark matter: relic density and direct detection jj] . 

In spite of its success and popularity, the MSSM with /2-parity still has two major prob- 
ems: the /i-problem jsf and the potential proton decay problem due to dimension 5 operators 
gI. The MSSM does not explain why the only mass parameter in the superpotential (/i pa- 
rameter) is at the electroweak scale instead of the fundamental (Planck or GUT) scale. 
Furthermore, i?-parity allows for dimension 5 lepton and baryon number violating operators 
such as QQQL and U'^U'^D'^E'^ which may still violate the bound on the proton lifetime. A 
natural solution to these problems would probably require that the MSSM be extended by a 
new mechanism or a new symmetry. The f/(l)'-extended MSSM (UMSSM) is a straight- 
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forward extension of the MSSM with a non-anomalous TeV scale Abelian gauge symmetry. 
It can solve the /i-problem and the dimension 5 operator problem very naturally with an ap- 
propriate charge assignment. An extra Abelian gauge symmetry is also suggested by many 
new physics scenarios including grand unification extra dimension j^, superstring 10 1, 
little Higgs theories [11], strong dynamics jl2|, and Stueckelberg mechanism |13 |. 

The UMSSM retains most attractive features of the MSSM, and also brings new aspects 
in relation to the dark matter problem. First, it extends the set of possible LSP dark 
matter candidates: the superpartners of the Higgs singlet S and the Z' gauge boson are now 



additional components of the neutralinos 



14| . The (predominantly right-handed) sneutrino 



also becomes a thermal CDM candidate if it couples to Z' ISj. Furthermore, a new gauge 



symmetry also naturally introduces exotic particles required to cancel gauge anomalies 
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When a neutral massive field is among the exotics, it could be a CDM candidate as well: 
it could be the LSP and stable due to i?-parity, but it is also "plausible that it is not the 

nn 

LSP, yet still is stable due to a new discrete gauge symmetry [17|, |18| , which may emerge 
naturally as a residual symmetry after the f/(l)' gauge symmetry is broken. Because of its 
gauge origin, the discrete gauge symmetry does not suffer from domain wall problem [l^ 
or symmetry violation by gravity effects pCl]. With the presence of two discrete symmetries 
very naturally introduced, the dark matter sector is enriched with coexisting multiple dark 
matter particles. Most previous studies on SUSY dark matter assume that the dark matter 
is made up of a single component. Then it is usually required that the LSP makes up 100% 
of the CDM in the Universe or, if the nature of the LSP is such that its annihilation rate is 
too high, SUSY dark matter is simply treated as a subdominant component of the CDM, 
without discussing the dominant component or providing the quantitative consistency with 
experimental constraints. 

In this paper, we consider the f/(l)'-extended MSSM with a residual discrete gauge 
symmetry {Z2 for simplicity) arising from the f/(l)' gauge symmetry. We consider a simple 
case where a SM singlet exotic is a stable particle under the new Z2 symmetry. Together 
with the LSP which is stable under /2-parity, there can be up to three dark matters due to the 
kinematics (see Appendix E]) • We also illustrate explicit examples for multiple dark matters 
that satisfy the relic density and direct detection constraints forming a viable texture of 
dark matters in the Universe. 

In Section [TTl we describe the possible remnant discrete symmetry of the f/(l)' extension 
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of the MSSM. In Section IIIH we discuss the rehc density and annihilation channels. In 
Section HVt we discuss the direct detection. In Section |Vl we perform numerical analysis 
on the relic density and the direct detection. In Section IVIt we summarize our results. In 
Appendix lAl we list the conditions where the multiple particles are stable. In Appendix [Bl 
we list the relevant lagrangians in our model. 

II. ;7-PARITY AND THE LIGHTEST [/-PARITY PARTICLE 

There are already at least two SM singlet chiral superfields in our model, the Higgs singlet 
S that spontaneously breaks the new symmetry with its vacuum expectation value (vev) 
and the right-handed neutrino A^^ that explains the neutrino mass. Higgs is not suitable for 
the odd particle under a parity since its replacement with a vev would destroy the discrete 
symmetry. Right-handed neutrino that forms LH2N'^ is not good either for the odd particle 
since lepton doublet also should have the odd parity and it is the active light neutrino that 
would be the lightest odd particle under new parity. We will consider a new SM singlet (X) 
- neither S nor A^'^ - that may be required from the anomaly cancellation condition, as our 
new dark matter candidate^. 

It is natural to think about possible remnant discrete symmetry jl?! . [isl from our 11(1)' 
gauge symmetry. A discrete gauge symmetry Zn would emerge if the discrete charges and 
the U{iy charges satisfy the following relation: 

z[S]=N, z[F,] = q[F,] + n,N (1) 

where z[Fi\ and q[Fi] stand for the f/(l)' charge and the Zn charge for a field Fi, respectively. 
z[Fi\, q[Fi\, Hi and A^ are all integers after an appropriate normalization of charges. Higgs 
singlet S is supposed to have q[S] = (mod A^) to keep the discrete symmetry after it 
is replaced by its vev (e.g. both SXX and (S) XX should be singlet under the discrete 
symmetry). 

We assume a Z2 {N = 2) discrete parity, which we shall call [/-parity (Up), as a remnant 
discrete symmetry of the f/(l)' gauge symmetry. We assign even Up to the SM fields, and 
odd Up to the X field. We assume either a scalar {(px) or a fermion (ipx) component 

^ For an example of a SM singlet required for the anomaly cancellation in a C/(l) '-extended MSSM, see 
Appendix III of Ref. [21|]. 
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(whichever hghter) of the superfield X is the Lightest tZ-parity Particle (LUP). The discrete 
charges are then 

q[X] = 1, g[SM fields] = (mod 2) . (2) 

The Z2 symmetry allows terms that contains X only in even numbers such as SXX. To 
satisfy the condition ([T]) of the discrete gauge symmetry, we need to assign f/(l)' charges as 
z[S] multiplied by an integer to the SM fields, but not to the X field. We assume SXX is 
the only term with X in the superpotential (thus, ipx is a Majorana particle), and we can 
assign either even or odd i?-parity (Rp) to (px (and the other to ipx)- It is obvious that the 
MSSM and the UMSSM superpotentials can allow such a charge assignment. 

To be an anomaly-free theory, the f/(l)' charges should satisfy the gauge anomaly condi- 
tions when whole particle spectrum is considered. We do not restrict ourselves to a specific 
particle spectrum, and just require the LUP to be a SM singlet exotic, while various exotic 
particles may exist with odd or even ^/-parity. We assume all exotics are vector-like for the 
SM gauge group so that the SM gauge group anomaly conditions are not altered. If a Zn 



is a discrete gauge symmetry originated from a U 
should satisfy the discrete anomaly conditions 
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y gauge symmetry, the discrete charges 
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where i runs through all SM fields and the exotic fields, and yi stand for the hypercharges 
normalized to integers for all fields. If there is only one exotic field X, the discrete charges 
of eq. ([2]) satisfy theses discrete anomaly conditions automatically except for the [gravity] ^- 
U{iy and [U{iy], which means we need more exotics with nonzero discrete charges. We 
will not give further constraints from the full gauge anomaly-free conditions keeping the 
possibility of various exotic fields open. 



Although SH2H1 term (i.e. effective fi term) greatly motivates the low energy scale of 
the U{iy resolving the /i-problem, it is technically involved to perform numerical analysis 
for all channels when the Higgs doublets and the Higgs singlet are coupled (through Yukawa 
term as well as D-teim). Therefore, for the sake of simplicity, we will limit our numerical 
analysis to the case of z[Hi] = z[H2] = with which the /i-problem is not solved in the 
way of Ref. [?!, i.e. by forbidding the original term H2H1 and allowing the effective n 
term SH2H1. The analysis of the UMSSM and the LUP dark matter in a more general 
setup where the Higgs doublets have non- vanishing f/(l)' charges and /i-problem is resolved 
will be discussed in other publication. In the limit of vanishing charges for Higgs doublets, 
SH2H1 term as well as SS*HiH* terms from the D-term will be forbidden making the Higgs 
singlet S an isolated physical eigenstate and Higgs doublets Hi, H2 the same as those of 
the MSSM. (For general properties of Higgses in the UMSSM in comparison with other 



models, see Ref. [2^.) The Z and Z' mixing will be zero, and the Z'-'mo {Z') and singlino 
(5*) components of the neutralino sector are also completely decoupled from the MSSM 
neutralinos since the mixing terms are proportional to z[Hi]. 

The superpotential is then given by, assuming three right-handed neutrino Dirac mass 
terms, 

W = HH2H1 + yuH2QU' + VdHiQD' + yNH2LN' + veHiLE' + \sXX . (9) 

Here, the Yukawa terms give the following relations among family universal^ charges: 

z\Q\ = -z[U'] = -z[D'] = riQzlS] (10) 
z[L] = -z[N'] = -z[E'] = nLz[S] (11) 
z[X] = -h[S] (12) 

where uq and are integers. With these charge assignments, Majorana neutrino mass 
term N'^N'^ (unless = 0) or SN'^N'^ are not allowed in general. We will assume neutrinos 
are Dirac throughout this paper^. Therefore, we have three free parameters (including z[S] 
not normalized to A^) for the f/(l)' charges in our numerical study: Uq, til, z[S]. 



^ Family non- universal U{iy charges may induce dangerous flavor-changing neutral currents [23|. However, 
the flavor-changing Z' may explain the discrepancies in rare B decays with appropriate parameter choices 
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^ Additional relativistic degrees of freedom would contribute to the ''He abundance, but be diluted by large 
Z' mass. It may explain the discrepancy of the ^He measurement 25|, |26 |. 



The dangerous dimension 5 proton decay operators QQQL, U'^U'^D'^E'^, and U^D'^D'^N'^ 
have commonly total f/(l)' charge {3nQ + nL)z[S] (up to overall sign) with the above charge 
assignment, and are forbidden if condition 

SuQ + nL^O (13) 

is satisfied. 

In the i?-parity violating supersymmetric model or any other model which does not have a 
stable dark matter candidate, the LUP can provide a good CDM candidate with an addition 
of the U{iy gauge symmetry. Depending on the detail of charges and spectrum, the U{iy 



may also ensure longevity of the proton without i?-parity (for example, see Ref. |21|). 



III. RELIC DENSITY 

The rehc density of dark matter is precisely measured as ^cdm^^ = 0.1111o;q}5 (2cr al- 
lowed range) by the WMAP+SDSS Present day relic density of a dark matter component 
is given by 

where Mpi is Plank mass, s{To) is the entropy density at present time, h is the normalized 
Hubble constant and the relic abundance Y{T) is defined as the number density divided by 
the entropy density. The abundance of the dark matter Y{T) can be calculated by solving 
the Boltzmann equation 



dT V 45 ^pi^^^)^^(TY-Y,^{TY) (15) 

where is an effective number of relativistic degree of freedom and Y^qiT) is the thermal 
equilibrium abundance, {av) is the thermally averaged annihilation cross-section times 
relative velocity. 

For our numerical calculation, we use micrOMEGAs [27] which performs rapid relic density 
calculation for 2 — 2 processes. We implement, on top of the MSSM fields and interactions, 
additional fields and interactions of the UMSSM as well as three right-handed Dirac neutrinos 
with negligible masses. The relevant lagrangians of the model are listed in Appendix [Bl 

We take vanishing soft trilinear terms limit {A = 0) in our numerical analysis. A non- 
vanishing AsxxS4>x(px soft term would separate masses of two components of complex 0x 



field, after S is replaced by its vev, wliicli we want to avoid for numerical simplicity. Witfi 
assumption of tlie zero cliarges for Higgs doublets (as discussed in Section [IT]), and sufficient 
mass splittings (so that co-annihilations are irrelevant), the annihilation channels for ipx 
and (px are given as follows. 



1. 


ipxi^x ■ 


ff {Z' mediated s-channel) 


2. 


i'xi'x ■ 


//* mediated s-channel, Z' mediated s-channel) 


3. 


^x^x ■ 


SS, Z'Z' (S mediated s-channel, ipx mediated t-channel) 


4. 


ipxi^x - 


SZ' [Z' mediated s-channel, ipx mediated t-channel) 


5. 


'4'xi^x ■ 


— > SS {Z' mediated s-channel, (px mediated t-channel) 


6. 


V'xV'x - 


Z'Z' {(j)x mediated t-channel) 


7. 




— s> SZ' {S mediated s-channel, (px mediated t-channel) 


8. 


i^xi^x ■ 


— > 4>x<Px ^' mediated s-channel, S and Z' mediated t-channel) 


9. 


^x^x ■ 


— ^ 4'x<t'x{.<P*x<P*x) {S and Z' mediated t-channel) 


10. 




-> // {Z' mediated s-channel) 


11. 




//* {S mediated s-channel, Z' mediated s-channel, 4 point interaction) 


12. 




SS {S mediated s-channel, (px mediated t-channel, 4 point interaction) 


13. 




-s> Z'Z' {S mediated s-channel, (px mediated t-channel, 4 point interaction) 


14. 




SZ' {Z' mediated s-channel, (px mediated t-channel) 


15. 




-> SS {Z' mediated s-channel, ipx mediated t-channel) 


16. 




Z'Z' {ipx mediated t-channel) 


17. 


(t^xct^x - 


SZ' (5 mediated s-channel, ipx mediated t-channel) 


18. 




'ipxi^x {S and Z' mediated s-channel, S and Z' mediated t-channel) 


19. 


4>x4>x{< 


Px^*x) ~^ i^xi^x {S and Z' mediated t-channel) 
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20. (px4>x{(px^*x) ^ 'S'5', Z'Z', SZ' {ipx mediated t-channel) 



The ibv annihilation includes features of Ref. 



Ref. 



and (px annihilation includes features of 



15| . Additional neutralinos Z' and S mix with each other, and we will call the mass 



eigenstates of them as x' from now on. 



IV. DIRECT DETECTION 



There are many direct detection experiments that attempt to detect the dark matter 
particles via nuclear recoil as the Earth passes through dark matter halo of our galaxy (for 



examples, see Refs. 
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30 



31 



32 



331]). Most experimental limits of the direct detection 
are given in terms of the cross-section per nucleon on the assumption of sole dark matter. 
Here, we will first discuss the general formalism for multiple dark matters. 

Let us introduce a parameter which parametrizes the fraction of the energy density of 
a dark matter i in our local dark matter halo and also in the whole universe. 

halo Q . 

^^ = 4ir-^ (16) 

PCDM "CDM 

where ^cdm^^ = Tli ^ih'^, and J2i = 1- The event rate per unit detector mass per unit 
time is given, from a dimensional analysis, by 

Rr^na {v) /Mn (17) 

where Mjq is a nucleus mass and is a number density of the dark matter. A precise 
differential rate [3i], for each component z, is given by 

^ = /"f; 2 F\En)UEn) (18) 

where Uo is the elastic scattering cross-section of a dark matter off a nucleus, fi^ = 
MM]\f /{M + Mtv) is the effective mass of a dark matter and a nucleus, and Vq ~ 220 km/s 
is the circular speed of the Sun around the galactic center (neglecting the relative velocity 
of the Sun and the Earth). 

Let us assume a simple form factor of F{Eji) ^ 1, and a pure Maxwellian speed distri- 
bution which gives 

Ti{ER) = expi-vl^Jvl) (19) 
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where f j,min = \J^^^ ^^e minimum velocity of a CDM that can produce an energy 
deposit above the detector threshold energy Et- Then the direct detection rate for each 
component is 

where Ei^ = 2fj,%VQ/M]\f, and the total detection rate is a sum of the contributions from each 
component. This total rate is the quantity constrained by the experiments. 

i? = J] i?, = C J] ^ exp i-ET/Ei,) < i?e.p (21) 

i i ' 

where C = 2p^^l,vo/{^/l^MN). 

Because of different masses of the CDM species, the usual constraint on the cross-section 
is not applicable in general for multiple dark matter cases. However, in some special cases 
where the relevant mass is effectively only one, those constraints can still be used. For 
example, let us consider a special case of where only one dark matter component (^) interacts 
with a nucleus while the others do not. The experimental bound on the detection rate 
R < i?exp can be rewritten for the cross-section 

eao < M exp {Et / Eo)C-^ R,,p (22) 

where e, (7o, Eq, M are those of the dark matter ^. It is the cross-section per nucleon that 
is usually compared to the experimental data, which is obtained by 

where = MMn/{M -\- Mn). The effective cross-section for the single detectable dark 
matter case is then given by cr„^efr — ^^ni ^-nd the experimental constraint is imposed on e(7„ 
instead of u^. 

Returning to our model, the dominant channels for the direct detection of ■(/'x and 0x 
are Z' mediated t-channels. If we had SH2H1 term, there would have been other channels 
such as Hi mediated i-channel (where is a general mixture of Higgs doublets and a Higgs 
singlet). For the MSSM LSP, we consider the hghtest neutralino (x*^)- 

The effective low-energy lagrangian for elastic scattering of the CDM and a quark is given 

by 
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where 



^^x-q = G z'z[X]zA[q]ipxl^lb'^xqit,lbq (25) 
£^,_, = iGz'z[X] {<P\d,<Px - {zvWiy + ZA[qWl,q) (26) 

£^o_, = + rf,3?7/.75X°g7^75g (27) 

and Gz' = gl,/Ml,, zv[q] = j{z[qL] + z[qR]), ZA[q] = K^^M - z[qL]). The fg term resuhs 
from Higgs exchange and squark exchange. The dg term is due to Z exchange and squark 
exchange. In the non-relativistic hmit the ipx interaction and the dg term are spin- dependent 
interactions, while the other interactions are spin-independent and dominate due to the 
coherence effect. 

With a choice of the MSSM-hke hghtest neutrahno (i.e. mass eigenstates of 5* and Z' 
are heavier than the hghtest mass eigenstate of B, W^, Hi, H2) and heavy squarks, the 
effective scalar coupling of a the neutralino dark matter to up-type and down-type quarks 
are approximately given by 

Hi=h,H ^« Hi=h,H 

Tfi^o^o = (A^4o cos a + N^o sin a) {N20 — tan 6^^^ A^io) (29) 
Thx^^o = sin a — N30 cos a) {N20 — tan 9wNiq) (30) 



_ g2mu cos a _ ^s^^sma 

2Mw sm p 2Mw cos p 

_ g2mu sin a _ g2md cos a 

■J- Huu— . n J-Hdd— n K'^^J 

2Mw sm p 2Mw cos p 

where a is the mixing angle that diagonalizes the CP even Higgs mass matrix, and tan (3 = 

m/iHi). 

In the non-relativistic limit, the spin- independent cross-section of a CDM and target 
nucleus is given by 

= dL^ZX, + {A-Z)Xn)' (33) 

TT 



with 
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where 



^x: Ap(^x) = 0, A„(^x) = (34) 
<Px ■■ Ap(0x) = Gz'z[X]{2zv[u] + zv[d]), A„(0x) = Gz'z[X]{zv[u] + 2zv[d]) (35) 

X° : Ap(x°) = 2/, = 2mp I frA + :^/tg ^ ) , A„(x°) = 2/„ ^ 2f, (36) 

with non-neghgible nucleoli parameters given by fr, = 0.12, /t^g = 0.84 (s-quark, heavy 
quarks contribution for both proton and neutron) 3J]. The spin-independent cross-section 
for ipx dark matter is zero. Then, with all three dark matters present, the experimental 
bound can be written as 

57 SI 

exp i-Er/Et) + ^ exp {-Er/Ef) < C-^R^^,. (37) 

V. NUMERICAL ANALYSIS 

In this section, we present the numerical results of the relic densities and direct detection 
cross-sections. For definiteness of our numerical analysis, we make following choices for our 
parameter values. 

1. We assume the U{1)' gauge coupling constant and charge assignments as following: 

9z' =9i = \ j\9Y (38) 

ng = 1 nz. = 5 z[S] = 0.3 (39) 

which correspond to 

^[gL]=^M=0.3 z[iL] = z[in] = 1.5 z[X] = -0.l5. (40) 

The above choice of ng and ul satisfies the condition (fT3|) to avoid dangerous dimen- 
sion 5 proton decay operators, and the relatively small f/(l)' coupling to quarks help 
avoiding the direct detection constraints for the (px dark matter. 

2. We choose the dark matters so that the lightest neutralino (x°) is always the LSP, 
and either the ipx (with Rp even) or the (px (with Rp even) is the LUP out of several 
possibihties. 
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3. We choose the parameter values as 

M2 = n < Ml (41) 

so that x° is wino or higgsino-hke. The wino or higgsino-hke LSP has been considered 
not a good cold dark matter candidate since their annihilation rate is too high to 
satisfy the measured CDM relic density. In the multiple dark matters scenario, they 
are actually preferable LSP dark matter candidates that can coexist with other dark 
matter. Since M2 and fi are the diagonal components of the chargino mass matrix, 
we usually get a light chargino (xf) also, and the co- annihilation between and Xi" 
becomes important for the relic density calculation of We need to scan M2 for 
analysis, but for other gaugino masses we choose 

Ml = 1000 GeV M3 = 1000 GeV Mi^ = 100 GeV (42) 

and choose Ma = 500 GeV, tan (3 = 40 and Ms = 800 GeV. 

4. We assume degenerate masses for squarks and sleptons and vanishing trilinear terms 

M,q„,,ks = 1100 GeV Msieptons = 600 GeV A,u = . (43) 

The large squark masses will make the squark contribution to the direct detection 
small. 

5. We assume all exotics (except X) that might be required for the anomaly cancellation 
are very heavy, and neglect their contributions to the relic density. 

In Figure [T], we illustrate the relic density and the effective spin-independent cross-section 
versus the LSP mass. We assume the LSP is x^ and the LUP is tpx with Rp even and 
M^^ = 440 GeV. 0x is assumed to be very heavy {M^^ = 1400 GeV) so that it is not 
stable. The green band is the 2a region measured by WMAP-I-SDSS. The relic density curve 
of the ipx is mostly flat since its annihilation cross-section is not sensitive to the M^o, and 
the CDM relic density is dominantly accounted for by the LUP dark matter for the given 
parameter values. The poles in the relic density curve of the x° are due to the resonances 
through (in the co-annihilation of and xf )? Z, h and H/A. The x° does not couple 
to the Z' and we flx Mz' = 1000 GeV for dark matter. M ± ~ M o and M o < 100 GeV 
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FIG. 1: (a) Relic density and (b) direct detection effective cross-section (e^ocj^o) versus the LSP 
dark matter (x^) mass in the presence of the LUP dark matter ^px■ We fix M^^ = 474 GeV and 
M(f)^ = 1400 GeV, Mz' = 1000 GeV. Dashed curve is the rehc density for tpx and sohd one is the 
total relic density of ipx and x^- Vertical line is the exclusion limit by the LEP bound on chargino 
mass. 



is ruled out by the LEP constraint on chargino mass of M^± > 104 GeV 35(] as indicated 
by the vertical line in the figure. The light Higgs mass also satisfies the LEP bound of 
Mh > 114 GeV over all range of the plot. We limit the range of the neutralino mass in the 
plot to be M^o < 500 GeV so that it is always the LSP where slepton masses are fixed at 
600 GeV. Although the relic density of the LSP dark matter varies with M^o, the total relic 
density is almost constant over the entire M^o. Therefore the finely measured relic density 
does not give a severe constraint on the LSP property except for its upper bound on the 
relic density. 

Since the ipx does not have a spin-independent cross-section with a nucleus (eq. ( !34l) ). the 
subdominant LSP dark matter is the only detectable dark matter by the spin-independent 
nuclear recoil experiment. The effective cross-section is proportional to the relic density, and 
the curve is valid for the entire M^o region, since the ^cdm^^ 7^ 0.1 region (M^o < 20 GeV) 
is already excluded by the LEP chargino mass bound. We assumed ^cdm^^ = 0.1 for 
calculation of = Qi/QcDM in the effective cross-section (e^oa^o) for the entire range of 
M^o. The green curves show the current CDMS and XENON experimental constraints on 
the direct detection. It is easy to satisfy the constraints since the flux of the detectable dark 
matter is smaller than that of the entire dark matters. 
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FIG. 2: (a) Relic density and (b) direct detection cross-section versus the LUP dark matter {(j)x) 
mass in the presence of the LSP dark matter x^- We fix M^^ = 1400 GeV, M^o = 150 GeV, 
Mz' = 500, 1000, 2000 GeV. Dashed curve is the rehc density for (px and sohd one is the total relic 
density of (j)x and x^- Thick black bands are where the total I^cdm/i^ satisfies the WMAP+SDSS 
limit. 

In Figure [21 we illustrate the relic density and the direct detection cross-section versus 
the LUP mass. We assume the LUP is (px with Rp even and the LSP is 4'x is assumed 
to be very heavy {M^^ = 1400 GeV) so that it is not stable. We show curves for Mz' = 500, 
1000, 2000 GeV to illustrate the dependence on Z' mass. The relic density curve of the x is 
fiat since its annihilation cross-section is not sensitive to the M^^. The CDM relic density 
is mostly dominated by the LUP dark matter. 

In the (px curves, we can observe effects of the annihilation channels we discussed in 
Section IIIII We will discuss several of them here. We see the Z' resonance poles at M^^ ~ 
Mz'/2 as expected. Among resonances through the S, only the (f)x4'*x ~^ x'x' for ^z' = 
500 GeV is open at M^^ ~ Ms/2 = 400 GeV. This channel is kinematically forbidden for 
Mz' = 1000 or 2000 GeV, since M^/ ~ Mz' ± Mi>/2 when Mz' > M^. This S resonance 
of Mz' = 500 GeV is the only region where (px relic density is smaller than that of x^- 
The new channel 4>x(p*x ~^ opens up for M^^ > Mgicptons = 600 GeV, making the relic 
density drop significantly. This dark matter annihilation channel where the final states 
are heavy supersymmetric particles (much heavier than the LSP) is a novel feature of this 
model that can not be found in the single LSP dark matter scenario like the MSSM. The 
effect of 4>x4>*x ~^ Z'^' is relatively small, but it is still noticeable for Mz' = 500 GeV at 
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M^^ ~ Mz'. The (i)x<P*x ^ SS channel for M^^ > Mg = 800 GeV is also distinguishable 
only for relatively light Z' case. This is because the dominant contribution in this channel 
comes from the F-term which is proportional to k"^ . The k is not an independent variable in 
our choice of input parameters, and it is small for large Z' mass. From eqs. fIBSP and flB7p . 
k is given as k = gz'z[S]M^^/Mz'- Due to the variety of channels, it is not difficult to find 
points with right relic density. Those points were marked as thick black band in the direct 
detection plot. 

Since the LUP is dominant {e^^ ^ 1), the effective direct detection curve is practically 
the cross-section of the LUP, except for the 5* resonance pole for Mz' = 500 GeV which does 
not satisfy the relic density constraint anyway. Therefore, we plot only the direct detection 
curve for the (px- Most part of the Mz' = 500 GeV curve with the choice of parameter 
values is excluded by the direct detection experiment, while the Mz' = 1000, 2000 GeV 
curves survive. The overall size of (px direct detection rate with TeV scale Z' is comparable 
to the current ongoing dark matter experiment. The cross-section decreases if Z' coupling 
to quarks gets smaller. If charges are quark-phobic (-2[q'L], ~ 0), a light Z' would also 
survive the current constraint by the direct detection and the Tevatron jsG] experiment. 

Besides the two examples we explored here, there are more possibilities such as all three 
dark matters coexist as well as only the LUP dark matter exists without the LSP. Also the 
LSP and LUP may form the same supersymmetric multiplet (dark matter supermultiplet) 
in a case, for example, ipx is the LSP and (px is the LUP. 

In the MSSM with a single LSP dark matter, the direct detection and the collider exper- 
iment is correlated and it is expected that both experiments detect the dark matter of the 
same property in the form of nuclear recoil and missing energy, respectively. With multiple 
dark matters, the expectations may change. 

For the neutralino LSP dark matter, the annihilation cross-section should be larger than 
the single dark matter scenario since it is responsible for only fraction of the CDM relic 
density. While the smaller flux will diminish the chance of the direct detection (in the 
case X is the only detectable dark matter), the increased coupling may enhance the chance 
of detection at the collider. The wino or higgsino-like LSP dark matter may also imply 
interesting phenomenology in the indirect dark matter search. For example, a large higgsino 
component could mean large capture rate in the Sun and correspondingly large neutrino 
fluxes probable in the IceCube experiment [37|. 
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For the LUP dark matter, the spin-independent direct detection rate by the LUP dark 
matter itself could be null {^jJx case) or large (0x case). The LHC/ILC may detect them 
if they are within the reach of colliders, ipx and (px couple only to Z' and S, and the 
prediction will change depending the Z' coupling to quarks and leptons. For example, the 
quark-phobic Z' will result in null direct detection rate even for (px dark matter and the 
LHC will not be able produce Z' resonance from the typical Drell-Yan process. 

Overall, with the presence of multiple dark matters with different sensitivity to the direct 
detection and to the colliders, the property of the dark matter measured in one experiment 
may not be consistent with that in the other experiment. 

VI. CONCLUSION 

TeV scale SUSY is a very well motivated new physics scenario. The minimal version 
of the supersymmetric SM (MSSM) may not be the correct reahzation of the TeV scale 
SUSY though, and its popular dark matter candidate (LSP) may not be the correct or 
full description of the dark matter. Various issues of the MSSM actually suggests it needs 
to be extended to include new ingredients such as new Abelian gauge symmetry f/(l)'. 
The concept of the supersymmetric dark matter may also be extended in the alternative 
supersymmetric models. 

We considered a residual discrete symmetry ([/-parity) which naturally emerges from the 
extension of the MSSM with the C/(l)'. It provides a new dark matter candidate (LUP) 
as a stand-alone alternative or coexisting complementary one to the usual LSP dark mat- 
ter. We showed that two well-motivated Z2 symmetries (i?-parity and [/-parity) can allow 
various interesting possibilities and numerically illustrated a few examples that satisfy the 
experimental constraints for the viable dark matter. 

The enriched dark matter properties suggest that the phenomenology for the supersym- 
metric dark matter (such as physics at the upcoming collider experiments) may be drastically 
different from the MSSM predictions. For example, the missing energy at the collider exper- 
iment may be originated from two or more massive stable particles, and it will be necessary 
to develop a technique that can distinguish that from one dark matter case. 
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APPENDIX A: STABLE PARTICLES 



There may be up to three stable massive particles when there are two parities. Here, we 
discuss it and categorize the possibilities in our model. With two symmetries [Hp and 
L/p), we can classify all component fields into four in terms of (-Rp, t/p). 

(+, +) : class A (— , +) : class B 

(+, — ) : class C (— , — ) : class D 

We call the lightest particle of each class as Aq, -Bq, Co, -Dq, respectively. (Aq is naturally the 
SM photon.) The stability of the lightest particles among i?p odd (LSP) and the lightest 
among Up odd (LUP) are guaranteed by Z2 symmetries. 

The minimal decay channels of Z-2, odd particles that conserve both parities are 

5o^Co + /^o, Co^fio + ^o, /^o^5o + Co. (Al) 

A non-minimal decay channel is, for instance, 

5o ^ oiC + 02D + n^A (A2) 

where Oj is an odd positive integer and nj is a non-negative integer. 
From the above, we can see that 

M^o < + M/Jq guarantees that Bq is stable. (A3) 

< M^g -|- M/Jq guarantees that Co is stable. (A4) 

< + guarantees that Dq is stable. (A5) 

If M^Q, M(7g, M/Jq satisfy a triangle relation (i.e. if all eqs. (]A3l IA41 lASP hold as Figure [3]), 
there will be three stable particles. Since C/p is originated from the gauge symmetry, which 
commutes with SUSY, any field F and its superpartner F share the same ^/-parity. The 
relevant fields are then classified as follows. 

1. quarks, leptons, gauge bosons, Higgses: class A (+, +) 

2. squarks, sleptons, gauginos, higgsinos: class B (— ,+) 

3. Hp even component of X: class C (+, — ) 
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FIG. 3: Triangle relations of masses (mi + m2 > for mi < m2 < m^) to have three stable 
particles under two Z2 parities. 

4. Rp odd component of X: class D (— , — ) 

There may be more exotics that belong to A and B, or C and D depending on their [/-parity, 
but for the sake of simplicity of the discussion we will consider only one exotic superfield 
X. Then, if triangle relation holds, both ipx and (px should be two components of the 
triangle (Co and Dq) and the lightest one among other superparticles (Bq) should be the 
third component of the triangle. The third component could be the LSP or the Next-to-LSP 
(NLSP) depending on the relative masses of Bq and Dq. 

To illustrate stability of three particles with triangle relation, let us assume = 
500 GeV, M^^ = 400 GeV and M^^ = 300 GeV, and assume x is the NLSP. We as- 
sign Rp even to (px and -Rp odd to ipx- i'x — ) is stable because it is the LSP; (px (+, — ) 
is stable because it is the LUP; x (~) +) is stable because its only decay channel x 4'x(px 
is closed kinematically. This corresponds to the Case III of Table [T] which contains the com- 
plete list of triangle relation cases. It is understood that when the triangle relation does not 
hold, stable massive particles will be less than three. 

APPENDIX B: ADDITIONAL LAGRANGIANS 

Here, we list the lagrangians from the f/(l)' symmetry and new field contents. We omit 
doublet Higgs and higgsino terms which has vanishing f/(l)' charges in our setup. //, and 
fa represent all the MSSM chiral fields including three Dirac neutrinos. 
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T 
1 


TT 
11 


TTT 
111 


IV 


V 


VI 


Rp odd X {Do) 


mi (LUP, LSP) 


mi (LUP, LSP) 


m2 (LSP) 


m2 (LUP) 


ma 


ma 


Rp even X (Cq) 


m2 


ms 


mi (LUP) 


ms 


mi (LUP) 


m2 (LUP) 


x{Bo) 


1713 


m2 


ma 


mi (LSP) 


m2 (LSP) 


mi (LSP) 



TABLE I: Possible triangle relations, x (^o) is the lightest among superparticles possibly except for 
the Hp odd component of X (Dq). A triangle relation mi +m2 > ma is assumed (mi < m2 < ma) 
among masses of il^x, (t>x and x- The LSP and the LUP can be one particle where the other two 
particles are NLSP and Next-to-LUP (NLUP). mi should be the mass of the LSP or the LUP (or 
both) , and ma can not be that of LSP or LUP. 



1. Fermion-Fermion-Z' 

A = -\9z'z[h]Z'^li,{l-i')f-\gz'z[fR]Z'^li,{l + i')f 

+\gz'z[X]Z'i'J^^^^'>'il^x + \gz'z[S]Z"''s^^^''S (Bl) 

2. Scalar-Scalar- 

A = -i9z'z[fL]Z'^{fldJr^-dJUr:)-igz'z[fR]Z'^^^^^ 

-igz'z[X]Z'^{cP\d^ct>x - d^rx<Px) (B2) 

3. Scalar-Scalar-Z'-Z' 

A = gl,z[h]'Z'^Z'Jlh + gl,z[fnrZ'^Z'J*Jn 

+gl,z[XfZ'^Z'^cj>\cl>x + glz[SfZ'^Z'^S*S (B3) 

4. Scalar-Fermion-Z' 

A = -^gz'z[fL]{z'{l-l'')m-rl{l-rl'>)Z'h) 
+^gz'z[fR] (J'{1 + f )//^ + 7(1 - j')Z'fn) 
-^gz'z[X] (J'{1 - 7')Vx0x + + l')Z'<l>x) 
-^gz'z[S] (J'{l-j')SS* + S{l + j')Z's) (B4) 
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5. D-term 



A = -^gh {zihUlh - z[fR]rjR + z[X]<f>*^<t>x + z[S]S*Sy (B5) 

6. F-term 

A = -k'S*S<l>*^<l>x - (B6) 

7. Yukawa 

(^(1 - 75)^x5 + ^(1 + 75)^x^*) (B7) 

8. Soft term 

Cs = (—MyZ'Z' - AsxxScPx^x + h.c\ - mlS*S - mU\ct>x (B8) 
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